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I.  IltJTRQDUQTXQN 

A.  Definitions 

1 .  The  Coagulation  of  particles  in  the  most  general  sense 
is  a  process  by  which  they  suffer  mutual  collisions,  which  may 
terminate  in  coalescence  (if  they  are  liquid  droplets)  or 
aggregate-formation  (if  they  are  solid). 

2.  The  dispersion  of  Particles  is  a  process  by  which  they 
occupy  increasingly  larger  portion  of  space  with  time. 

B.  Significance 

Both  particles  coagulation  and  dispersion  are  very 
significant  in  real  life. 

As  an  example,  when  disseminating  fibrous  or  platelet-like 
material  as  an  aerosol  in  the  field,  the  process  of  coagulation 
(which  mainly  occurs  in  the  early,  high-concentration  stage  of 
the  chain  cf  events)  will  determine  particle  size  and 
(frequently)  complicated  shape. 

Hence,  this  process  will  affect  light-scattering  properties  of 
the  cloud  and  its  deposition-rate  on  various  obstacles  in  its 
path;  it  will  also  influence  dispersivity  of  the  cloud. 

Aerosol  cloud  dispersion  will  determine  particle  concentration 
at  various  space-time  points  (x,t)  in  the  field. 
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Aerosol  coagulation  ensues  by  several  mechanisms  which 
operate  in  laminar  and  turbulant  flow  fields: 

i.  Brownian  coagulation  (in  laminar  and  turbulent  fields). 

ii.  Gravitational  coagulation  (in  laminar  and  turbulent 
fields),  which  occurs  due  to  the  different  fall  velocity  of  the 
particles . 

iii.  Gradient  coagulation  (in  laminar  and  turbulent  fields), 
which  occurs  due  to  the  passing  of  particles  over  each  other  as 
they  are  carried  by  their  stream-lines  and  thus  intercept  each 
other . 

iv.  Inertial  coagulation  (in  laminar  and  turbulent  fields) 
which  occurs  due  to  the  different  inertial  response  of 
particles  to  the  velocity  field  of  the  medium  (air),  thus 
causing  the  emergence  of  a  new  relative  velocity  between  the 
former. 

We  will  be  interested  in  processes  within  a  turbulent 
atmosphere  and  in  non-spherical  particle.^,  which  is  a 
completely  new  research-line. 
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2 .  Di.spersion 

Aerosol  dispersion  ensues  mainly  under  turbulent  conditions; 
it  is  caused  by  the  stochastic  response  of  particles  to  the 
(turbulent)  velocity  fluctuations  of  the  medium  (air)  and  thus 
is  a  major  property  of  the  turbulent  field. 

The  stochastic  response  depends  on  particle  size,  shape  and 
turbulence  characteristics. 

We  will  be  interested  in  the  stochastic  response  of  non- 
spherical  particles  to  a  turbulent  field,  which  has  been  only 
scarcely  studied,  since  coagulation  may  occur  only  when  the 
cloud  of  aerosols  is  concentrated. 

II.  PREVIOUS  STUDIES 

We  denote  by  m  ,  i=l,2,3  the  concentration  of  single, 
double,  triple  particles  etc. ,  formed  by  the  process  of 
coagulation . 

The  Collision  Coefficient  C^j  is  defined  by  the  relation-' 

Q  L;  = 

Q-^j  is  the  collision  rate  of  particles  i  and  j  with 
concentrations  n^.  and  nj  . 

The  concentrations  nk  satisfy  the  following  set  of 
equations : 


Jli  ^  y-7 


[2] 
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A .  Brownian  Coagulation 


Smoluchowski  (1)  in  his  pioneering  work  on  Brownian  Coagulation 
in  a  still  medium  assumed  that  Cij=C.  The  solution  to  the 
set  of  equations  [2]  under  this  assumption  is: 


h 


;  t)  =  . 

HCnoh 


[3] 


Where  no  is  the  initial  concentration  of  identical  particles 
and  n  is  the  total  number  of  particles.  Eq.  [2]  was  later 
solved  numerically  without  the  assumption  of  equal  coeficients. 


B.  Turbulent  Coagulation  of  Spherical  Particles 


Here  three  cases  have  been  treated  differently. 

1 .  Asymptotic  Case  A 

Spherical  particles  which  are  much  smaller  than  Kolmogorof f ' s 
length  scale 

If  -  kinematic  viscosity  of  medium  (air)  0.15  cm2 /sec 

e  -  energy  dessipation  rate  per  unit  mass  5-200  cm2 /sec^ 
in  the  lower  atmosphere. 

In  the  lower  atmosphere  =  0.1-0. 5  cm 


Saffman  &  Turner  (2)  obtained  a  formula  for  Ci  j  which  takes 
into  account  the  three  mechanisms  ;  the  inertial  response  of 
particles  to  the  flow  of  air,  gravitation  and  gradients  of  the 
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velocity  of  the  air. 

a.anf  Cd-z/j-/  (x,-x/  ( + 

Here  R,  :  =  r.  +  r^  ,  r,  being  the  radius  of  the  i-th  particles 

*•7  7  *> 

formed  by  i  original  particles. 

/  -  is  the  density  of  air 
_fp  -  density  of  particles 

T;  =  2  /p  /?>w  “  ^he  relaxation  time  of  a  Stokesian 

particle . 

M  -  dynamic  viscosity 

g  -  gravitational  acceleration 

-  mean  square  of  total  time  derivative  of  th€  air 

01/ 

velocity . 


This  last  quantity  may  be  expressed  in  terms  of  known 

quantities  of  the  flow  according  to  Batchelor  (3): 

— .  1/z. 


(puj*  =  ijftVi/) 


If  gravitational  and  inertial  effects  are  neglected,  the 
expression  obtained  is: 

Ci;,  =  i-J  (6/j/J 

(It  differs  by  a  small  factor  from  the  expression  obtained  from 
[5]  upon  substituting  and  due  to  approximations 

made  in  deriving  it).  If  one  assumes  that  =  2C  one  obtains 
for  the  densities  the  same  solutions  as  in  Brownian 
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coagulation,  with  the  characteristics  of  the  turbulent  field 
entering  into  the  constant  C. 

2.  Asymptotic  case  B 

Spherical  particles  which  are  much  larger  than  Kolmogorof f ‘ s 
length  scale.  Abrahamson  (4)  obtained  the  following  formula: 


v'  is  the  r.m.s.  velocity  of  the  particle.  It  is  related  to 
the  r.m.s.  velocity  of  the  fluid  u'  by 

Equation  [7]  is  valid  for 

Ir"  >6aviA^(  Wpj/e 

3 .  Intermadia^e  Case  C 

This  is  the  case  for  which  particle's  sizes  are  between  the 
sizes  dealt  with  in  cases  A  and  B  .  A  theory  for  spherical 
particles,  belonging  to  this  case  has  been  developed  by 
Williams  &  Crane  (5). 

III.  PRESENT  STUDY 

A.  Aiin.  Study 

Inspite  of  the  fact  that  most  real  aerosol  systems  are  composed 
of  non-spherical  particles,  very  little  has  been  done  so  far 
with  respect  to  the  investigation  of  simultanuous  coagulation 
and  dispersion  of  such  particles.  This  study  is  consecutive  to 
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those  of  Krushkal  &  Gallily  (6), (7)  on  particles  orientation 
which  were  supported  by  CRDC,  U.S.  Army. 

As  pointed  out  in  the  research  proposal  the  aim  of  this  study 
was  to  produce  a  method  which  would  be  able  to  predict 

a)  particle  shape  and  size  distribution  and 

b)  Particle  dispersion, 

both  under  the  influence  of  a  sufficiently  characterized 
turbulence.  This  study  was  limited  to  the  case  of  particles 
much  smaller  than  the  Kolmogoroff  size. 

Simultanuous  dispersion  and  coagulation  of,  initially 
monodispersed ,  elongated  ellipsoides  in  a  simultated  turbulent 
flow  have  been  calculated  up  to  the  stage  of  two  particles 
collision.  The  original  aim  of  the  study  was  to  calculate  also 
collisions  of  a  single  particle  with  an  aggregate  and  too  take 
account  of  the  hydrodynamic  interaction  between  particles.  Due 
to  the  death  of  the  chief  investigator,  Prof.  Isaiah  Gallily, 
these  goals  have  not  been  achieved. 

B.  General  Method 

The  coagulation-dispersion  problem  was  treated  by  a  Monte  Carlo 
technique  in  which  an  ensemble  of  non-spherical  particles  was 
subjected  to  the  effects  of  a  fluctuating  velocity  and 
fluctuating  velocity  gradients  in  a  turbulent  air  Due  to  the 
fluctuating  factors,  the  particles  suffered  a  fluctuating 
translation  and  rotation  which  caused  them  to  coagulate  and 
disperse . 

The  pafticles  were  assumed  to  be  uniformly-sized  prolate 
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ellipsoids  of  revolution  of  a  large  aspect  ratio.  They  were 
also  taken  to  be  much  smaller  than  the  typical  Kolmagoroff  eddy 
size  {  y  li  ) 


In  the  Monte  Carlo  calculations,  the  tragectories  of  the 
particles  were  calculated  from  their  equation  of  motion. 

The  inertial  memory  of  the  partii les  was  taken  into  asccount 
through  their  mass  and  moment  of  inertia 


TTi^y  =  + /X  k  ‘  (u-v) 


[9] 


=  L,  th. 


[10] 


where  K  and  L  on  the  r.h.s.  of  these  equation  are  the  Oberbeck 
and  Jeffery’s  solutions  given  in  (8),  (9)  and  (10). 


L,  -  |!>n/U»bc  r  (b  -C)/ +  (fe  -  U)J] 

3(bY.4c‘^;r 

u,  V  are  fluid's  and  particles  velocity  respectively,  a,b,c, 
are  the  semi-axes  of  the  ellipsoidal  particles. 


[11] 


[12] 


[13] 


b>a=c,  3=b/a 
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(3o 


a*b 


c 


Cb'+>U 


OO 

To'  — 

/  (.c'+7i)A 


f  is  a  strain  tensor  component 


f  = 


DZ  I 


"I  is  a  vorticity  component. 

1  ::  f  Oy] 

^  ~  32  ^ 

u,  V,  w,  f  and  ^  are  expressed  in  a  system  of  coordinates 
attached  to  the  ellipsoid  and  parallel  to  its  principal  axes. 

It  was  assumed  that  during  each  time  interval  of  particle 
motion,  the  velocity  gradient  was  kept  constant,  though  it 
changed  from  one  interval  to  another  according  to  a  probability 
function,  truncated  around  its  zero  mean. 


IV.  RESULTS 

A.  Turbulent  field  properties 

For  the  turbulent  field  we  took  the  following  values: 

e  =  103cm2/sec3,  V=  0.15cm2/sec. 


Thus  for  the  typical  Kolmogoroff  length  we  had: 

)^=  4.3«10-2cm 

The  Kolmogoroff  typical  velocity  was  3.5  cm/sec  and  the  typical 


time-1 . 2» 10~ 2  sec. 
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The  velocity  gradients  within  each  eddy  were  assumed  to  be 
normally  random  and  constant  during  a  time  interval  At;^  .  The 
velocity  ij  is  given  by 

Ulx  +  J*)=  UU)+  ti4: 

The  normal  distribution  of  the  velocity  gradients  was  truncated 
at  i  4  ^  ,3  given  by  the  following  relations  between  the 

gradients  and  typical  values  of  the  turbulent  field. 


e/15 


(i=j=k=l) 


-  V 


-e/30 


2e/15 


(i=k=j=l  or 
k=o=l=i) 

(k=l=i= j ) 

(all  other  combinations) 


[15] 


B.  Method  of  Calculations 


The  particles  were  given  randomly  distributed  initial 
translational  and  rotational  velocities  with  zero  mean.  Since 
the  relaxation  time  of  the  particles  was  much  shorter  than  a 
time  interval,  the  particles’  translational  velocity  was, 
practically,  the  velocity  of  the  air.  the  rotational  velocity 
was  calculated,  at  each  time  step,  from  equation  [13].  At  the 
end  of  each  time  interval  t=10“3sec.,  new  gradients  were 
generated  at  the  location  of  the  particle  according  to  [15]. 
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The  new  translational  velocity  of  the  particle  was  calculated 
according  to  [14]  and  the  new  rotational  velocity  according  to 
[13]  and  [10],  with  the  velocity  gradients,  appearing  in  [13], 
taken  at  the  new  point.  During  each  time  interval  collisions 
were  checked  at  intervals  of  10" < sec.  in  2-D  calculations  and 
2*10~5sec.  in  3-D  calculations. 


C.  2-D  Calculations 


Calculations  were  performed  first  in  two  dimensions  in  order  to 
check  the  method. 

The  particles  were  prolate  ellipsoids  of  revolution  with  semi 
axes  b  and  a,  chosen  according  to  a  log  normal  distribution. 

The  geometrical  mean  length  -  <2b>  =  2»10'3cm 

g.x.m.s.  =  1.2*<2b> 

The  geometrical  mean  width  -  <2a>  =  2»10''>cm 

g.r.m.s.  ^  =  1.2»<2a> 

V 

Truncation  was  at  ±  0.15  of  the  mean  size 
Particles  density  fp=l.g/cm3 
The  time  intervals  were 

(  lO'^^sec  for  collision  checking 


nt;  - 


10-3  sec 


for  constant  velocity  gradients 


Total  time  of  calculation  :  SAt^ =3* 10- 2  sec . 

Initial  inter-particle  distance  =  10-2cm. 
Initial  particle  concentration  =  10^ /cm3 
Initial  number  of  particles  =  100. 


Fig.  1  and  Fig.  2  show  two  typical  numerical  experiments.  In 
these  experiments  the  particles  were  separated  after  each 
collision  and  the  coagulation  continued.  In  each  of  the 
experiments  the  number  of  particles  was  64 .  In  the  first  one 
(Fig.  1)  27  collisions  were  counted,  in  the  second  -  12. 

D.  3-D  Calculations 


As  in  the  2-D  calculations  the  particles  were  prolate 
ellipsoids  of  revolution  of  log-normally  distributed  sizes.  We 
aimed  at  obtaining  relations  between  the  size  and  the  initial 
concentration  of  the  particles  and  the  number  of  collisions. 
Therefore  calculations  were  performed  with  two  sets  of 
particles’  sizes  and  concentrations. 


Sizes 

i . 
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regular  particles 

mean  length  -  <2b> 
g.r.m.  3. 

mean  width  -  <2a> 
g.  r . m. s . 
large  particles 

All  linear  dimensions  are  twice  as  large. 


4 • 10- 4  cm 
1.2-<2b> 

0 . 8« 10“ ^  cm 
1 . 2*<2a> 


Truncation  was  at  ±  0.15  of  the  mean  size. 


i.  regular  concentration 


initial  number  of  particles  N=  125 

inter  particle  initial  distance  X=  10“ 3 cm 

initial  number  concentration  no=103cm“3 

ii.  high  concentration 

initial  number  of  particles 
inter  particle  initial  distance 
initial  number  concentration 
In  addition  the  effect  of  the  initial  number  of  particles  was 
checked. 


N=  125 
X=  S^lO-'^cm 
no  =8* 103  cm“ 3 


Initial  .number 

i.  N=64 

ii.  N=125 


particles  density 


=lg/cm3 


time  intervals 


2« 10“ 5  sec  for  collision  checking 


At-,  ^ 


10“ 3  sec  for  constant  velocity  gradients 


Total  time  of  calculation  :  Zat,  =  2« 10“ 2  sec. 

1..— Pigp.ergipn 

A  typical  picture  of  the  dispersion  of  the  particles  is  shown 
in  Fig.  3.  It  shows  that  after  2-10“2sec  the  particles  are  so 
dispersed  that  it  is  unlikely  to  expect  many  more  collisions  to 
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occur . 


2.  Collisions 

The  number  of  collisions  depends  non-linarily  on  the  initial 
concentration  of  the  particles.  According  to  Saffman  &  Turner 
(2)  the  number  concentration  of  pairs,  as  a  function  of  time 


n^=  j  xr  5-.I- RMe/i/f't 


[16] 


R  is  the  radius  of  initial  spheres. 

For  nol<<l 

n2  2^  no^X  [17] 

In  our  calculations  ikjuf  -  1 . 2«  10“  2  sec  1  .  The  values  of 
X/t,  based  on  half  the  length  and  half  the  width  of  the 
paeticle  and  on  the  initial  number  concentration 
no =8« 109  cm" 2 ,  are: 


a.  For  R=4«10"5cm 

I/t=4-10"i 5 

[18] 

b.  For  R=2*10"'’cm 

t/t=5*10-i 3 

[19] 

In  Fig.  4  typical  results  of  numerical  experiments  are  shown 
with  particles  of  size  (i),  number  concentration 
no =8* 109  cm" 3  and  initial  number  N=125.  It  shows  that  most 
collisions  occur  in  the  first  3*  10" ■’sec  and  that  the  total 
number  of  collisions  after  this  time  is  30-40.  According  to 
[17]  and  [18]  the  time  needed  for  this  number  of  collisions  is 
-  10"  ■’sec  and  according  to  [19]  -  10"  8  sec.  The  longer  time 
needed  in  our  calculations  is  due  to  the  strong  dispersion 
which  was  not  accounted  for  in  (2). 
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In  Fig.  5  similar  results  are  shown  for  larger  particles  of 
size  (ii)  and  initial  concentration  no =109  cm' 3. 

3 .  Effect  of  Concentration 

In  Fig.  6  the  number  of  collisions  is  given  for  three  initial 
concentrations  (the  graph  shows  the  average  over  great  number 
of  numerical  experiments).  The  lowest  concentration  - 
no =1 . 25* 108 cm' 3  seems  to  exhibit  a  too  low  number  of 
collisions  for  any  conclusion  to  be  drawn  from  it.  The  average 
slope  ai ,  corresponding  to  initial  concentration 
no  =109  cm' 3  in  the  first  2*10~'<sec,  is  ai  -  5/2 
whereas  the  slope  a2 ,  corresponding  to  no=8*109cm-3  is 
a2  25/2  so  that  a2/ai=5. 

Similar  results  are  shown  in  Fig  7.  Here  the  initial  number  of 
particles  was  N=64.  The  ratio  of  the  slopes  of  the  two  graphs 
taken  at  t=l .  4"  10" sec  is  a2/ai=6.25.  It  should  be  noted 
that  the  leveling  off  in  the  number  of  collisions  for  N=125 
occurs  after  t=2»10"4sec  whereas  for  N=64  after 
t=l . 4* 10' 4  sec .  This  is  due  to  the  lower  number  of  particles. 
These  results  are  in  disagreement  with  [16]  and  [17],  from 
which  a  ratio  of  a2/ai=64  is  to  be  expected.  This 
disagreement  is  a  consequence  of  the  small  number  of  particles 
introduced  initially  in  the  experiment  and  of  the  rapid 
dispersion  shown  above.  Nevertheless,  the  strong  dependence  of 
the  number  of  collisions  on  the  initial  concentration  is  still 


obvious . 
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4.  Effect  of  Size 

Larger  particles  are  heavier  but  the  mass  does  not  enter  into 
the  expression  [16].  Also  it  does  not  enter  in  our  calculations 
since  the  relaxation  time  was  small  compared  to  the  time  step. 
Therefore  inertial  effects  were  not  significant  in  the  range  of 
sizes  we  have  applied.  Thus,  the  only  parameter  that  determines 
the  results,  if  the  initial  number  of  particles  N  is  large 
enough,  is  the  ratio  of  the  spacing  between  particles  and  their 
sizes.  An  increase  by  a  factor  of  2  of  each  linear  dimension  of 
a  particle  should  therefore  be  equivalent  to  an  increase  by  a 
factor  8  of  the  initial  concentration.  In  Fig.  8  the  number  of 
collisions  is  shown  for  regular  particles  of  size  (i)  and  for 
large  particles  of  size  (ii).  The  initial  number  is  N=125.  In 
Fig.  9  comparison  is  made  between  regular  particles  of  size  (i) 
with  initial  concentration  no=8»109cm-3  and  between 
larger  particles  (ii)  and  initial  concentration 
no=109cm"3.  The  similarity  between  the  results  assures  us 
of  the  validity  of  the  method  of  calculations. 


5.  Effect  of  Initial  Number 

Our  calculations  were  performed  with  a  small  number  of 
particles  due  to  computational  limitations.  The  increase  of  the 
initial  number,  i.e.  the  volume  of  equal  concentration  of 
interacting  particles,  results  in  an  increase  in  the  number  of 
collisions.  This  is  roughly  explained  by  the  fact  that  if  two 
equal  volumes  of  particles  are  put  together,  particles  from  the 
two  volumes  will  interact  through  the  common  surface.  Thus,  the 
number  of  collisions  will  be  greater  than  the  sum  of  collisions 
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in  the  two  separated  volumes.  The  effect  of  increase  in  N  is 
shown  in  Fig.  10.  It  should  be  noted  that  not  only  the  number 
of  collisions  increased  with  the  increase  of  N,  but  also  the 
time  interval  at  which  high  rate  of  collisions  occur  is  larger. 
Thus,  with  N=64  the  time  at  which  the  graph  of  collisions 
levels  off  was  t^=^  1. 4*  10"  ^  sec  whereas  for  N=125 
t— 1 . 4» 10~  4  sec . 


A  numerical  model  was  constructed  for  the  simulation  of 
dispersion  and  coagulation  of  a  cloud  of  non-spherical  aerosol 
particles  in  turbulent  conditions.  Coagulation  was  observed  in 
the  early  stages  of  the  dispersion  of  the  cloud.  It  appears 
from  the  calculations  that,  in  a  small  cloud,  dispersion  is  the 
dominant  process  from  the  early  stages  of  the  formation  of  the 
cloud  and  coagulation  stops  after  a  short  time.  Nevertheless, 
at  the  time  coagulation  occured  it  exhibited  a  behaviour 
similar  to  what  is  expected  from  existing  theories  but  modified 
by  the  strong  dispersion  of  the  cloud.  The  intention  of  the 
chief  investigator,  Prof.  I.Gallily,  was  to  perform  numerical 
simulations  with  a  much  larger  number  of  particles,  to  examine 
the  behaviour  of  the  system  with  particles  of  different  aspect 
ratios  and  to  take  account  of  the  particle-particle 
hydrodynamic  interaction.  Due  to  the  death  of  Prof.  Gallily 
these  additional  goals  were  not  achieved. 
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Fig.  4  -  Number  of  collisions 
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NP^initial  number,  RO=initial  number  concentration 


125  PARTICLES  -  CONCENTRATION  EFFECT 


T*E-4  sec 

Number  of  collisions  for  3  initial  concentrations  RO 
particle  size  (i),  initial  number  NP=125 


1  ,  E'”4  sec 


Fig.  10  -  Number  of  collisions 

NP  -  initial  number  of  particles 

HC -  -  concentration  (ii),  particle  size  (i). 


